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Abstract 

The many-body state of carriers confined in a quantum dot is controlled by the bal- 
ance between their kinetic energy and their Coulomb correlation. In coupled quan- 
tum dots, both can be tuned by varying the inter-dot tunneling and interactions. 
Using a theoretical approach based on the diagonalization of the exact Hamiltonian, 
we show that transitions between different quantum phases can be induced through 
inter-dot coupling both for a system of few electrons (or holes) and for aggregates 
of electrons and holes. We discuss their manifestations in addition energy spectra 
(accessible through capacitance or transport experiments) and optical spectra. 
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1 Introduction 

Semiconductor quantum dots (QDs) are a formidable laboratory for next- 
generation devices and for the actual realization of some key Gedankenexperi- 
mente in many-body physics [1-3]. Indeed, the number of electrons and holes 
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in the QD can be controlled very accurately, and almost all relevant parame- 
ters influencing their strongly correlated states, like confinement potential and 
coupling with magnetic field and light, can be tailored in the experiments. 
The additional possibility of tuning the coupling between QDs enriches their 
physics and the possible applications. 

From the point of view of fundamental physics such coupling extends the 
analogy between quantum dots ("artificial atoms" [4]) and natural atoms, to 
artificial and natural molecules. The tunability of coupling among QDs allows 
to explore all regimes between non-interacting dots and their merging into a 
single QD; many of those regimes are precluded to molecular physics. 

One of the peculiarities of QDs with respect to other solid state structures 
consists in the partial decoupling of a few degrees of freedom from all the 
others, which is due to the discrete nature of the spectrum [1-3]. The actual 
exploiting of such a feature largely depends on the capability of integrating 
arrays of QDs, thus increasing the number of degrees of freedom that one 
can address with precision and coherent manipulation. This is precisely the 
strategy pursued by the semiconductor-based solid state implementations of 
quantum computation [5]. 

In general and basic terms, the tuning of inter-dot tunneling allows to modify 
the relative position of the single-particle levels, thus inducing phase transi- 
tions in the many-body ground states and different degrees of spatial corre- 
lations among carriers. Manifestations of these phenomena in systems formed 
by carriers of only one type, whose ground and excited state properties are ac- 
cessible through addition energy spectra, have been predicted. Here we point 
out that similar effects are expected to occur also for systems formed by both 
electrons and holes. We also show that, in spite of the obvious differences, 
strong similarities appear in the analysis of electrons and electron-hole sys- 
tems, and a unified theoretical description is in order. Basically, a competition 
emerges between two trends. On one side an atomic Aufbau logic, where carri- 
ers tend to occupy the lowest single-particle states available, thus minimizing 
the kinetic energy and the total spin, at the (energetic) cost of reducing spa- 
tial correlation among carriers. At the opposite extreme we find an enhanced 
degree of spatial correlation among carriers, which occurs through the occu- 
pation of orbitals other than the lowest. This implies an enhancement of the 
kinetic energy and a reduction of the repulsive one, and results in electron dis- 
tributions maximizing the total spin (Hund's rule). The balance between these 
two trends depends on the spacings of the single-particle levels involved, and 
these are precisely what can be settled by controlling the inter-dot tunneling. 
When carriers of opposite charge, different effective masses and tunneling pa- 
rameters come into play, the competition between both trends becomes even 
more delicate. 
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Predictions of the actual ground and excited states of the many-body sys- 
tem thus require a careful theoretical treatment including all carrier-carrier 
interactions. Since the number of carriers in the dot can be controlled and 
kept relatively small, we can proceed through direct diagonalization of the 
exact many-body Hamiltonian, with no need to make a priori assumptions on 
the interactions. On the contrary, the results are a useful benchmark for the 
validity of the most common approximations for these systems. 

We find that different quantum phases correspond to different regimes of inter- 
dot coupling both for a system of few electrons (or holes) and for aggregates of 
electrons and holes, with various possible spatial configurations and the forma- 
tion of different possible "subsystems" of inter-correlated particles. Besides, 
due to the negligible electron-hole exchange interaction in heterostructures 
such as GaAs, the two kinds of carriers can be treated as distinguishable par- 
ticles. Therefore spatial correlation among electrons and holes does not arise 
from the Fermi statistics: it needs instead the entanglement between the or- 
bital degrees of freedom associated to holes and electrons, and turns out to 
depend only indirectly on the spin quantum numbers S e and SV 

After a brief summary of the state of the art in theoretical and experimen- 
tal work on coupled dots (Sect. 2), in the following we describe the general 
Hamiltonian and solution scheme (Sect. 3). We then come to the results for 
electron- (Sect. 4) and electron-hole systems (Sect. 5). The trends leading to 
different quantum phases are discussed in detail, together with their nature in 
terms of spin and spatial correlation functions. 



2 Experimental and theoretical background 

Early experimental and theoretical studies focused on electrostatically-coupled 
dots with negligible inter-dot tunneling [6] . Here we consider artificial molecules 
[7], where carriers tunnel at appreciable rates between dots, and the wavefunc- 
tion extends across the entire system. 

The formation of a miniband structure in a one-dimensional array of tunnel- 
coupled dots was demonstrated more than a decade ago [8]. After that, the 
first studies considered "planar" coupled dots defined by electrodes in a two di- 
mensional electron gas. In these devices the typical charging energy was much 
larger than the average inter- level spacing, hence linear [9] and non-linear [10] 
Single Electron Tunneling Spectroscopy (SETS), obtained by transport mea- 
surements at different values of the inter-dot conductance, could be explained 
by model theories based on capacitance parameterizations [11]. Early studies 
also considered simple model Hamiltonians (usually Hubbard-like) with ma- 
trix elements treated as parameters [12]. Blick and coworkers clearly showed 
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the occurrence of coherent molecular states across the entire two-dot setup, 
analyzing transport data [13] and the response to a coherent wave interfer- 
ometer [14]. The tuning of coherent states was also probed by microwave 
excitations [15], and coupling with environment acoustic phonons was studied 
[16]. Planar coupled dots were also used to cool electron degrees of freedom 
[17], to measure the magnetization as a function of the magnetic field [18], and 
to study the phenomenon of "bunching" of addition energies in large quan- 
tum dots [19]. The so-called "vertical" experimental geometry was introduced 
later: it consists of a cylindrical mesa incorporating a triple barrier structure 
that defines two dots. Sofar, evidence of single-particle coherent states in a 
AlAs/GaAs heterostructure has been reported [20], while in AlGaAs/InGaAs 
structures clear SETS spectra of few-particle states have been observed as a 
function of the magnetic field B and of the inter-dot barrier thickness [21]. 

A relevant part of theoretical research has addressed the study of few-particle 
states in vertical geometries, within the framework of the envelope function 
approximation. The two-electron problem was solved, by means of exact di- 
agonalization, in different geometries by Bryant [22] and by Oh et al. [23]. 
Systems with a number of electrons N > 2 at B m in cylindrical geometry 
have been studied by several methods: Hartree-Fock [24], exact diagonaliza- 
tion for N < 5 [25], numerical solution of a generalized Hubbard Model for 
N < 6 [26] and for N > 12 with a "core" approximation suitable for the 
weak-coupling regime only [27], density functional theory [28]. 

Palacios and Hawrylak [29] studied the energy spectrum in strong magnetic 
field and negligible inter-dot tunneling with various methods (N < 6), and es- 
tablished a connection between the correlated ground states of the double-dot 
system and those of Fractional Quantum Hall Effect systems in double layers. 
In this perspective, Hu et al. [30] studied collective modes in mean-field theory, 
Imamura et al. [31] exactly diagonalized the full Hamiltonian at strong B and 
different values of tunneling (N < 4), Mart in- Moreno et al. [32] considered the 
occurrence of canted ground states. Also the far-infrared response of many- 
electron states was analyzed with various techniques [33] . Another interesting 
issue is the relation between quantum and "classical" ground states [34] as the 
radius of the dot is enlarged, when electrons arrange to minimize electrostatic 
repulsion because the kinetic energy is quenched [35,36]. 

The electronic properties of planar dots have also been studied theoretically, 
through a variety of techniques: configuration interaction or analytical meth- 
ods with various approximations [37,38], or density functional theory for larger 
values of N [39]. The infrared [40] and the thermoelectric [41] response were 
considered. Systems of coupled QDs are also among the most promising can- 
didates for the implementation of semiconductor-based quantum information 
processing devices: some of the current proposals identify the qubits with ei- 
ther the spin [42] or the orbital degrees of freedom associated to the conduction 
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band electrons in QDs. 



Research on few-electron systems in double dots is thus a new field in very 
rapid growth, with increasing focus on the possible quantum phases and how 
they can be driven by artificially controllable parameters such as inter-dot 
coupling, magnetic field, dot dimension. The study of such phases is expected 
to add insight into the physics of double layers, e.g. the conditions for Wigner 
crystallization, and of strongly correlated systems in general. The amount of 
experimental data on many-body states in artificial molecules is still limited, 
but the whole bunch of spectroscopic tools currently available (linear and non- 
linear transport, Raman spectroscopy) is now beginning to be employed (see 
e.g. SETS spectra in the B-N space, Ref. [43] in this issue) and should allow the 
direct verification of theoretical predictions and a more general understanding 
of the basic phenomena and trends. 

Also the optical properties of coupled QDs depend both on the confinement 
of electrons and holes and on the effects of correlation among these carriers. 
In spite of their importance, however, such correlation effects are still largely 
unknown. From the experimental point of view, cleaved-edge overgrown sam- 
ples have been used [44] , but self-organized quantum dots are most commonly 
employed for optics. Their stacking was demonstrated [45], and the splitting 
of the excitonic ground state in a single artificial molecule was studied as a 
function of the inter-dot distance. The lines in the photoluminescence spec- 
tra were explained in terms of transitions among excitonic states obtained 
by single-particle filling of delocalized bonding and anti-bonding electron and 
hole states [44,46]. When a few photoexcited particles are present, however, the 
correlations induced by the carrier-carrier Coulomb interactions play a crucial 
role [47] ; single-particle tunneling and kinetic energies are also affected by the 
different energetic spacings of electron and hole single-particle states. The cor- 
related ground and excited states will thus be governed by the competition of 
these effects, not included in previous theoretical descriptions of photoexcited 
artificial molecules [50]. A detailed understanding of exciton and multiexciton 
states in coupled semiconductor QDs, however, is of great interest for the de- 
velopment of the optical implementations of quantum-information processing 
schemes, starting from the identification of well characterized qubits [51]. The 
possibility of complete optical control over the computational space formed 
by interacting excitons in quantum dots has recently been demonstrated in 
Refs. [52] and [53]. We therefore expect that a systematic investigation of 
trends in the many-body phases of coupled dots will be actively pursued for 
systems of few electrons, and extended to systems of photoexcited electrons 
and holes in the near future. 
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3 Many-body states of iV electrons and holes 



In the following we focus on the motion of few electrons and holes confined in 
two coupled quantum dots. Our primary interest is in the correlated nature 
of ground and excited states of the interacting system. 

Hereafter we consider a simplified model where, within the envelope function 
and effective mass approximations, two coupled identical vertical dots are 
described by a separable confining potential V(g,z) = V(g) + V(z), with 
V(g) = ^m*u 2 g 2 an in-plane parabolic potential [g— (x, y), m* is the effective 
electron (hole) mass, ujq the characteristic frequency] and V(z) a double square 
quantum well along the z direction (see Fig. 1). Each well (of width L and 
barrier potential height Vb) corresponds to a dot; the coupling between the 
two dots is controlled either by varying the inter-dot distance d (width of the 
inter-dot barrier) or the height of the inter-dot potential barrier. To vary d 
implies to consider differently grown devices. The full many-body Hamiltonian 
7i (in zero magnetic field) is the sum of the single-particle terms H^°\r) = 
—h 2 V 2 /(2m*) + V(g, z) and of the two-body Coulomb interaction terms: 
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Here n r is the dielectric constant of the semiconductor medium, and the sub- 
script e (h) refers to electrons (holes). Effective masses, characteristic frequen- 
cies, and details of the double well entering differ for electrons and holes. 

We choose this geometry for two reasons: firstly, experimental devices whose 
behavior can be described by this model are currently studied by several 
groups, allowing for precise tailoring of the dot geometry, strong spatial con- 
finement, and hence observation of spectral features beyond the simple Coulomb 
Blockade behavior (e.g. in SETS spectra). Secondly, the cylindrical vertical 
geometry, contrary to in-plane devices, has the richest degree of symmetry, 
which is particularly helpful to theoretical work both in reducing the size of 
Hilbert space sectors and in analyzing electronic configurations. Specifically, 
7i is invariant under any rotation in the spin space (the total spin S and 
its projection S z are therefore conserved), rotation around the z axis in real 
space (conservation of the ^-component of the orbital angular momentum M) , 
inversion with respect to the geometrical center of the system (parity conser- 
vation). In complete analogy with Molecular Physics [54] and for each species 
of carriers we introduce a spectroscopic notation to classify electronic terms, 
namely eigenstates of Ti: 2S+1 M 9tU . Here g (u) stands for even (odd) parity 
and M takes the labels S, n, A, . . . standing for M — 0, 1, 2, . . . Actually, a £ 
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term is also invariant under reflection with respect to a plane passing through 
the symmetry axis: in this case the notation takes the form 25+1 where ± 
labels the sign change under reflection [55]. 

We are interested here in the evolution of the ground and excited states as the 
inter-dot distance d is varied. This feature shows a remarkable difference be- 
tween artificial and natural molecules: in the latter the inter-nuclear distance is 
almost fixed, controlled by the nature of bonding, while in the former it can be 
tuned by adjusting electrods or by growing different sample devices. Ground 
and excited states can be probed by several kinds of spectroscopies. Theoret- 
ically, once the energy spectrum is known after numerical diagonalization of 
H, it is quite easy to compute the relevant observable quantities. 

A considerable achievement has been obtained by transport spectroscopies, 
like single-electron capacitance tunneling spectroscopy [56] or SETS [57] for 
the ground state, or non-linear tunneling spectroscopy [58] for the excited 
states. In a transport experiment the chemical potential fJ,(N) of the double- 
dot is measured as the number of electrons N is varied charging the system 
one electron by one. In fact, from the experimental value fi(N) one can infer 
information about the ground state, being ji (N) = E (N) — E (N — 1), with 
E (N) the ground-state energy of the iV-body system [59]. Our theoretical 
strategy is straightforward: we compute the ground state energies E (N) at 
different values of N, and from these the chemical potential \i to be compared 
with the spectra. Single-dot far-infrared spectroscopies [1-3] are unsuitable to 
probe the relative motion of electrons and hence their correlation, because light 
only couples to the center-of-mass motion (generalized Kohn theorem) [60]. 
This is also true for a system of vertically coupled quantum dots with cylin- 
drical symmetry, as long as the in-plane confinement potential (orthogonal to 
the symmetry axis, e.g. the growth direction) is parabolic and the polariza- 
tion of light is in the same plane. However, this limitation does not hold for 
two-photon processes like Raman scattering, where density fluctuations can 
excite collective modes of the interacting system [61]. 

Finally, optical spectroscopy allows the study of few-particle states including 
electrons and holes. In the lowest order the light-semiconductor coupling is 
associated either to the absorption of a photon and to the promotion of an 
electron from the valence to the conduction band or to the reversed process, 
which is accounted for by a Hamiltonian of the form — E • P, where E is 
the electric field and P the material polarisation [62]. Within the rotating- 
wave and dipole approximations the luminescence spectrum for a QD initially 
prepared in state |A) can be computed according to Fermi's golden rule: 

L a {u) oc ]T | (P CT )v, A | 2 5(E X + hu- E y ); (2) 

A' 

here (-Po-)a,a' are the dipole matrix elements corresponding to the transition 
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between states A and A' (through removal of one electron-hole pair) and the 
the creation of a photon with helicity a — ±. 



4 Few-electron system 

In this section we study the system of interacting carriers of the same species, 
e.g. electrons. Let us start from the simplest case, that is the two-electron 
molecule. A theorem due to Wigner [63] guarantees that the ground state is 
always a singlet if time- reversal symmetry is preserved [64] : however, dramatic 
alterations of the energy spectrum and wavefunction are driven by the inter- 
dot distance d and the characteristic dot radius £ = (h/m*^^) 1 ^ 2 . 

This is shown in Fig. 1: in panels (a) and (b) we plot the total ground state 
kinetic (E k ) and Coulomb (V ee ) energy [65], respectively, vs d. Coulomb in- 
teraction mixes up different configurations (i.e. Slater determinants) which 
contribute with different weight to the ground state. Besides (Ek) and (V ee ) 
for the true few-particle groundstate \ip) (diamond symbol), in Fig. 1 we also 
show the corresponding data of three prototypical states [66]: |1) = \a g |, a g [) 
(singlet), |2) = (\<j g ],a u [) - \a g [,a u \))/\/2 (singlet), |3) = (\a g ],a u j 
) + \ a gli cr ut))/v^ (triplet). Note that the difference between the (identical) 
kinetic energies of states \2) and |3) and that of |1) amounts exactly to the 
energy splitting A sas between the single-particle states a u and a g [Fig. 1 (a)]. 
This quantity decreases exponentially as d increases and as the probability of 
the tunneling through the potential barrier goes to zero. While singlet and 
triplet states |2) and |3) have identical kinetic energy, the latter state is ener- 
getically favored as the interaction energy is concerned. The splitting in (V ee ) 
between |2) and |3) appearing in Fig. 1 (b) is an exchange energy, namely the 
consequence of the antisymmetry of the total wavefunction for particle per- 
mutations. The behavior of the ground state \ip) partly resembles that of the 
state |1), but shows significant deviations due to the mixing of configurations. 

The arrangement of electrons is naturally visualized by computing density 
functions in real space. However, both the single-particle density and the usual 
radial pair correlation function g(g) plotted in the xy plane depend only on 
the relative distance, due to the cylindrical symmetry of the system. Hence, 
we follow Ref. [67] and calculate the "angular" spin-resolved pair correlation 
function 

P s , so (q,z; Qo,zo) = A SjS0 - Q)$( z{i) - z)5 s wj(g u) - £o)<K^ (i) - z )S s 

(3) 



where (...) denotes the expectation value on a given state, the subscript s 



refers to spin, and A StSo is a normalization factor, such that / dgdz dgo d^o Ps,s 
(g, z; go, zq) = 1. One electron with spin sq is fixed at the position (go, zq), while 
the other at (g, z) with spin s is varied: thus P s ,s (q, z; go, z ) is proportional to 
the conditional probability of finding the second electron given that the first 
one is fixed. This allows for observation of the relative spatial arrangement of 
electrons and of the angular correlation. The spin-independent quantity is the 
total pair correlation function P(g, z; go, zq), normalized as 



P(g, z; Qo, Zq) = 



]T N S (N S - l)P SyS0 + ]T N s N So P SyS0 



/N(N-l), (4) 



being N s the number of spin-s electrons. 



In Fig. 1 (c) we plot the function p{z) = JjdxdyP (x, y, z; x , y , z ) , showing 
how the fixed position of one electron (represented by the black circle) affects 
the spatial distribution of the other one along the symmetry axis z, for an 
inter-dot distance of 1 nm. The state |1) clearly exhibits no spatial correlation 
among the two carriers: the placing of one electron in one quantum dot (QD) 
does not change the probability of finding the other one in any of the dots. 
In the case of the singlet state |2) the spatial distribution of one electron 
is peaked around the other (fixed) one: the two particles tend to occupy the 
same QD. Opposite trends apply to the triplet state |3). Again, the true ground 
state shows a mixed character: p{z) has its biggest peak in the "unoccupied" 
QD, but there is a finite probability for the double occupancy on the same 
QD. The average values of the Coulomb energy (V ee ) [Fig. 1 (b)] clearly reflect 
such behaviors. The curve referring to |1) slowly decreases, because the value of 
{V ee ) corresponding to two particles in different QDs diminishes as d increases. 
A fortiori, (V ee ) decreases for the triplet state: the electrons are always in 
different QDs. The Coulomb energy is less affected by the inter-dot distance 
in the case of the state |2), because (V ee ) is mainly due to intra-dot interaction 
(both carriers in the same QD): the slight increase of (V ee ) depends on the 
growing localisation of the particles within a QD. 

The terms contributing to the Hamiltonian H of Eq. (1) scale differently with 
the characteristic length of the confinement potential i§. the kinetic one goes 
like ~ £q 2 , while the interaction one like ~ Eq 1 . For small dots, the kinetic 
term dominates and the system is Fermi-liquid like: here the ground state is 
determined by the successive filling of the empty lowest-energy single-particle 
levels, according to Aufbau atomic theory. As £ increases, the electrons become 
more and more correlated and arrange to minimize Coulomb repulsion, up to 
the limit of complete spatial localization (reminiscent of Wigner crystallization 
in 2D). Even if for N = 2 the ground state is always a 1 S 9 term as £q is 
varied, nevertheless we can gain further insight into the correlation dynamics 
by analyzing the two-body wavefunction. 
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In the inset of Fig. 2 we plot the total pair correlation function P(g, z; go, zq) 
for the N = 2 triplet state. Here, go an d z o are se t equal to the average value 
of the in-plane radius and the maximum along z of the single-particle density, 
respectively. In addition, z is fixed at the position of the second, symmetric 
maximum of density in the symmetry-axis direction: the resulting contour plot 
is the value of P in the xy plane (in units of £q). This represents the probability 
of finding one electron in the plane of one dot, given that the second electron 
is fixed on the other dot. The other plots in Fig. 2 show P(g, <p, z; go, Po, z o) 
vs the azimuthal angle ip: all other parameters g, g , z, z are fixed, with 
g — g - When p> = <p , the position coincides with that of the fixed electron, 
and the probability P has a minimum (zero in the triplet case with z = zq) 
due to the Pauli exclusion principle. As tp is varied, the position follows a 
trajectory like the thick circle in the inset, starting from the bullet locating 
the other fixed electron in the xy plane. After a 27r-rotation, we are back 
in the starting point. These plots are a kind of "snapshot" of the angular 
correlation, as we freeze the motion of one electron. Figure 2 is organized 
in two columns, corresponding to the singlet ground state and to the triplet 
excited state, respectively, for different values of ftw (d — 1 nm). Solid lines 
refers to the case z = z , namely electrons on the same dot, while dashed lines 
to z 7^ z , i.e. electrons on different dots. When fku is very large (40 meV, top 
row), the curves P(g,z; go,zo) are almost flat. This flatness implies that the 
motion of the two electrons is substantially uncorrelated, except the effect of 
Fermi statistics. In fact, in the triplet case, the probability of measuring two 
electrons on the same dot is negligible, and this holds at any value of hw . On 
the contrary, in the singlet state there is a finite probability of measuring two 
electrons on the same dot in the ground state. As hujo is reduced (20 meV, 
middle row), angular correlation is turned on. This can be seen by the increase 
of the peak- valley ratio in the angular correlation function. The position of the 
maximum corresponds to 7r, i.e. the two electrons repelling each other tend to 
be separated as much as possible. This trend is even clearer at small values of 
hu (3.5 meV, bottom row). Plots of the second row (Hljo = 20 meV) should 
be compared with the corresponding plot of Fig. 1 (c): the first ones illustrate 
how electrons correlate in the xy plane, the second one how they arrange along 
z. 

Note that in the first column of Fig. 2 as hu is decreased the probability 
of measuring two electrons on the same dot increases up to the limit when 
it equals the probability of measuring electrons on different dots {fkoo = 3.5 
meV). This is due to the different ratios between two fundamental energy 
scales: the harmonic oscillator inter-level separation ftwo and the energy differ- 
ence A sas between antisymmetric and symmetric double-well wavef unctions: 
if hwo ^> A sas , the inter-dot tunneling is negligible with respect to the kinetic 
energy of the intra-dot motion, and the dots are almost quantum mechani- 
cally decoupled. In the opposite limit, the system is coherent, and it makes 
no difference between measuring one electron on one dot or on the other one, 
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since the system behaves as a unique dot, doubled in size. 

Let us now turn to discuss the case N > 2. We choose a particular set of 
parameters, namely m* = 0.067m e , k t = 12.4, L — 12 nm, V = 250 meV, 
huo = 5.78iV -1//4 , corresponding to a set of experimental devices currently 
under study [21]. The parameterization of Hujo(N) is meant to mimic the 
effect of the gate voltage on the electrostatic confinement potential V{q) [28]. 
We exactly diagonalize the Hamiltonian 7i of Eq. (1) for N < 6, using up to 
32 single-particle orbitals. The convergence is checked controlling a cutoff on 
the average energy of the Slater determinants entering the computation. Our 
code uses the ARPACK package [68] and isolates Hilbert space sectors with 
S and S z fixed, contrary to usual Lanczos approaches. 

Figure 3 shows the calculated ground state energy vs d for 3 < N < 6. As d is 
varied, one or two transitions between ground states of different symmetries 
occur. Specifically, while there is only one transition between two different 
electronic terms for N = 3, two transitions take place for N = 4 and iV = 5. 
The intermediate phase for N = 4 exists only in a very narrow range of d (~ 
0.01 nm) in the neighborhood of d — 3.45 nm. For N = 6, again only two 
phases exist. However, at the intersection point of the 1 S 9 and 3 £ 9 terms, the 
excited state 3 Tl g is almost degenerate in energy. 

These transitions can be understood analyzing the many-body wavefunction 
of the different ground states [26]. In the bottom panel of Fig. 3 we focus on 
the N = 4 case and we schematically depict the major-weight Slater deter- 
minant corresponding to each phase. The key point is that, as d is decreased 
from the value of 4.5 nm, the "energy-gap" A sas between "bonding" and "anti- 
bonding" orbitals (i.e. symmetric and antisymmetric solutions of the double 
well along z) changes, from the limit of decoupled dots [labeled as c) in fig- 
ure] to the strong-coupling limit [labeled as a)]. In the c) case, the first-shell 
molecular orbitals 0a g and 0a u are almost degenerate and well separated in en- 
ergy with respect to the second shell, hence they are filled with four electrons 
giving the configuration Oo^Ocr^, i.e. two isolated dots with the first orbital 
shell completely filled. In the opposite limit, at small values of d [a) case], 
the bonding mini-band made of 0a g , 07r+, 0n~ single-particle orbitals is much 
lower in energy that the anti-bonding one. The ambiguity of how to fill the 
lowest-energy orbitals, due to the degeneracy of 0-7T+ and 07r~ levels, is solved 
consistently with Hund's first rule [57], i.e. the two open-shell electrons oc- 
cupy each orbital with parallel spin (the configuration being OagOn^On'), in 
such a way that exchange interaction prohibits electrons from getting close, 
minimizing Coulomb repulsion. This configuration is characteristic of a single 
dot, doubled in size [57]. In the intermediate phase b), the antibonding 0a u 
level is almost degenerate with the bonding Ott^ and 07r~ levels: while the 
first two electrons occupy the lowest-energy orbital 0<r 9 , the remaining two 
arrange again to maximize spin, consistently to a "generalized" Hund's first 
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rule. However, now there are three levels almost degenerate, and we find that 
the ground state configuration is 0<7 2 0cr u 07r+: according to Hund's second rule, 
also the total orbital angular momentum M is maximized, to minimize the 
interaction energy (the higher m, the smaller the Coulomb matrix element 
between single-particle levels). A similar reasoning applies to the transitions 
at iV ^ 4. 

Let us now focus on the N = 5 case in Fig. 3. As d increases, the ground 
state sequence is 2 I1 M — * 4 S M — > 2 I1 U , that is the 2 Tl u term appears twice, 
corresponding to a continuous energy curve that crosses twice the 4 S n term. 
However, if we examine which Slater determinants mainly contribute to 2 H U , 
we find that the relevant configuration at small d (I = 0<r 2 07r+ 2 07r~) differs 
from that at large d (II = 0<T 2 0o" 2 07r+). Moreover, the slope of the curve in 
the two regions is different, mainly due to the change in the balance between 
bonding and anti-bonding levels occupied: 5:0 for I, 3:2 for II, which controls 
the dependence of the overall kinetic energy on d (see the previous discussion) . 
This change in the "character" of the 2 U U term (i.e. the ratio between the 
weights of configurations I and II) is found to be continuous with d. We plot 
also the first excited state for the 2 n M symmetry (dashed line in the N = 5 
panel): clearly this curve anti-crosses the 2 n„ ground state. Analyzing the 
slope and character of this excited state in the small- and large-rf regions, we 
find an inverted behaviour with respect to the ground state: now the relevant 
configuration at small d is II while that at large d is I. The overall behavior can 
be understood as a consequence of the Wigner-von Neumann theorem, i.e. that 
intersection of terms of identical symmetry is forbidden [69]. Therefore, the two 
2 n„ terms anti-cross, while 2 Ii u and 4 S M terms can freely cross and bring about 
ground state transitions, belonging to different irreducible representations of 
the symmetry group of Ji. An analogous anti-crossing between ground and 
excited state for the symmetry is depicted in the N = 4 panel (solid and 
dashed line, respectively). 

Results of Fig. 3 should be compared with those obtained by means of exact 
diagonalization of a generalized Hubbard model [26], by density functional 
theory [28], and by Hartree-Fock method [24]. In all these works the window 
in rf-space at which the ground state 3 n 9 at N = 4 occurs is much larger and 
a ghost additional intermediate phase at N = 6 (corresponding to the excited 
state 3 n 9 in our calculation) appears. Therefore our results, that agree well 
with data obtained up to N < 5 by exact diagonalization in Ref. [25], clearly 
demonstrate the importance of correlation beyond mean-field approaches. The 
interacting electronic system is so correlated in regimes of realistic parame- 
ters of the devices that it is very difficult to obtain quantitatively reliable 
results with any approach but configuration interaction. This point was al- 
ready stressed, for single quantum dots, in Refs. [70] and [71]. 

To further characterize different ground states vs d, in Fig. 4 we plot the 
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spin-resolved pair correlation function P-[,-\(q, z; Qo, zq) for N = 4 for values 
of d corresponding to the three phases previously discussed. The position of 
a spin-up electron (qo,z ) is fixed in one dot as in Fig. 2, and the contour 
plots of the top (bottom) row, with z = z (z ^ z ) fixed, correspond to the 
probability of measuring another spin-up electron on the same (other) dot in 
the xy plane. The right column refers to the x E g term at d — 3.9 nm: there 
is only one "free" spin-up electron, and we can see that the probability of 
measuring it on the same dot as the fixed electron is negligible, while the 
probability distribution on the other dot depends only slightly on the position 
of the first fixed electron. Therefore the two dots are quantum mechanically 
decoupled, each one filled with two electrons in the lowest shell. The motion 
of electrons in the two dots is almost uncorrelated in the xy plane. In the 
opposite limit of small d (left column, d — 3.1 nm, three spin- up electrons), 
the coupling is so strong that it makes no difference either measuring the 
electron on one dot or on the other, i.e. the contour plots are identical, and 
the system forms a coherent, strongly bound molecule. The fixed electron is 
"dressed" by its exchange and correlation hole, i.e. it repels other electrons 
that are at small distances. The middle column {d = 3.5 nm, three spin-up 
electrons) shows that in the intermediate phase 3 U g the dots are coupled with 
a weak degree of coherence, namely the probabilities of measuring the electron 
on the two dots are sensitively different. Planar correlation in the same dot 
is important, while it is negligible for motion on different dots. The above 
classification of phases is also suitable to iV ^ 4. 

In this section we have shown results for the electron energy spectrum up 
to N < 6. It is straightforward now to compute the SETS linear spectrum, 
and the comparison with very recent experimental data [43] show remarkable 
agreement in many respects. Results for B > will be presented elsewhere 
[72]. 



5 Few electron-hole pairs 

We next consider systems of interacting carriers composed of an equal number 
of electrons and holes. Let us start by considering a single electron-hole pair 
(exciton) and the way in which its ground state depends on the width of 
the barrier. As the inter-dot distance d increases the splitting between the 
energies of the bonding (a g ) and anti-bonding (a u ) states decreases both for 
electrons and for holes. The energetic cost associated to the promotion of the 
two particles from the bonding to the anti-bonding states becomes smaller 
and comparable to the gain in Coulomb energy arising from the correlation 
of their spatial distributions along z. In Fig. 5 we plot the functions p{z) of 
electrons and holes, associated to given positions of the other carrier, at inter- 
dot distances of d — 1 nm (a) and d = 3 nm (b); the two insets represent 
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the contributions to the electron-hole ground states of 1 1) = \a e g f\ T) 
and |2) = \a e u T)- The decreasing tunneling goes with an increasing 

correlation (electron more localized around the hole and vice versa) and an 
increasing contribution from the |2) state. The slight difference between the 
plots associated to electrons and holes at each inter-dot distance are due to 
the differences in the barrier heights (400 meV for electrons and 215 meV 
for holes) and in the effective masses (m* = 0.067 mo and m* h = 0.38 mo) of 
the two carriers: as a consequence electrons tunnel more than holes and tend 
to be less localized within one QD. It is worth noticing that at d — 3 nm 
the electronic tunneling still induces a pronounced splitting between the two 
delocalized bonding and anti-bonding single-particle states (e& = 35.11 meV 
and e a = 37.52 meV). In spite of this, due to Coulomb correlation and to the 
reduced tunneling of holes, the energetic value and the spatial distribution 
of the excitonic ground state closely resembles that of an exciton in a single 
QD; besides, the splitting between the ground and the first excited states 
is negligible [73]. In other words the "excitonic tunneling" is suppressed at 
smaller inter-dot distances than the electronic one. 

If the double dot is occupied by two electrons and two holes, both attractive 
and repulsive interactions are present. Intuitively one would expect carriers 
with identical charge to avoid each other and carriers of opposite charge to 
look for each other: the interplay between such trends is directed by the values 
of d and by those of S e and SV In Fig. 6 we compare such correlations for 
two different values of the electron and hole spin quantum numbers and for 
d — 1 nm. The singlet-singlet lowest state (S e — 0, S h — 0) is characterized 
by a small correlation between the electrons, (a), and by a more pronounced 
one between holes, (c). Analyzing the eigenfunction associated to this state, 
one observes that approximately only the electron single-particle state a e g is 
(twice) occupied, while for holes strong contributions arise from both bonding 
and anti-bonding states. As already mentioned, this difference between the 
behaviors of the two carriers depends on the fact that a gain in Coulomb en- 
ergy has a greater kinetic cost for electrons than for holes. Besides, the spatial 
distribution of holes, (e), is not affected by the position of the electrons [74]. 
Let us now compare these correlation functions with the corresponding ones 
associated to the triplet-singlet configuration (S e — 1, S h — 0). Again the cor- 
relation among electrons and holes is negligible, (f): the two electrons and the 
two holes subsystems can thus be understood to a good extent independently 
from each other. As in the case of the prototypical state |3), that the state of 
the two electrons here resembles, the probability of finding two electrons in 
the same QD is negligible, (b). Such spatial separation of the electrons induces 
a more pronounced separation for the holes too, as shown by the flattening 
of the smallest peak in the left well [see Fig. 6 (b), (d)]. The differences be- 
tween the two spin configurations are even more dramatic at bigger inter-dot 
distances. In Fig. 7 we show the same correlation functions at d = 3 nm. 
The triplet-singlet configuration [Figs. 6 (b), (d), (f)] shows the same features 
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already observed at d — 1 nm; here the two electrons (holes) are perfectly lo- 
calized in different QDs, due to the suppressed tunneling. The singlet-singlet 
configuration [Figs. 6 (a), (c), (e)] instead has undergone a transition to a 
phase in which all carriers are localized in either QD (due to symmetry). If 
the position of one of the four particles is fixed in one QD, all the others are 
localized in the same one. This somehow surprising effect can be explained in 
the following way: in a mean field picture, due to the substantial similarity 
of the electron and hole wavef unctions, the localization of the two excitons 
in two different QDs or in the same one makes no difference with respect to 
the Coulomb energy because there is cancellation of terms of opposite sign. 
If correlation comes into play, however, the localization of all particles in the 
same QD gives rise to the so-called "biexcitonic binding energy" AE (which 
is defined as the difference between twice the energy of the excitonic ground 
state and that of the biexcitonic one). Specifically, the binding energy AE is 
due to the correlations among the x and y directions: as in the case of the two 
electrons in Fig. 2, such correlations become strongly effective and lower the 
Coulomb energy when particles are localized in the same QD. 

The comparison between the correlation functions corresponding to equal 
spin configurations at different inter-dot distances shows a trend similar to 
that of the two electrons alone. The population of the anti-bonding states in- 
creases with decreasing bonding-antibonding splitting, thus allowing a more 
pronounced spatial correlation between identical carriers. Such dependence is 
particularly clear for the singlet states, while for the triplet ones a high de- 
gree of correlation is already guaranteed by the permutational symmetry of 
the few-particle wavefunction (i.e., by the fermionic nature of electrons and 
holes) . 

The described behaviours are reflected in the values of the different contribu- 
tions to the mean Coulomb energies of each spin configuration (Fig. 8). Let us 
start by considering the three spin arrangements (S e , Sh) = (1, 1), (0, 1), (1, 0). 
The contributions to the Coulomb energy associated to the electron-electron 
interaction (a) all monotonically decrease with increasing d: the two electrons, 
each in a different QD, get more and more far apart. If S e — 1 the spatial 
separation of the two electrons is a direct consequence of the permutational 
symmetry of the wavefunction; if S e — the same effect arises from a proper 
linear combination of different Slater determinants and from the correspond- 
ing occupation of the electronic o e u orbitals (which in turn depends on the 
tunneling: from which the different slopes of the S e — and S e — 1 curves). 
Analogous behaviours are seen to occur with respect to the hole-hole interac- 
tion (c). The main difference as compared to the previous case is the higher 
degree of correlation among such carriers in the Sh = case, already at small 
inter-dot distances. The trends of the electron-hole Coulomb (b) interactions 
are hardly distinguishable one from the other. The monotonic decrease (in 
modulus) of (V e h) reflects that of the interaction energy among carriers lo- 
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calized in different QDs when they get far apart. The plots associated to the 
singlet-singlet configuration (continuous lines), however, show a transition to- 
wards a phase in which all carriers are localized in a same QD, already put in 
evidence in Fig. 7. The absolute values of all Coulomb terms correspondingly 
go through an abrupt increase for values of d in the range between 2 and 2.5 
nm. Let us finally consider the total energies (d). The lowest singlet-singlet 
state turns out to be the system's ground state at any inter-dot distance. For 
d < 2 nm the S e = and S e — 1 configurations are degenerate with respect to 
the value of S h ; the difference between the total energies follows that between 
the two-electrons triplet and singlet states. As d increases, the energies of all 
spin configurations but S e — Sh — asymptotically tend to a value which is 
twice the energy of the excitonic ground state in a single QD. The energy of 
the singlet-singlet state instead tends to that of the biexcitonic ground state in 
a single QD: the difference between these two asymptotic values is the already 
mentioned biexcitonic binding energy AE. 



6 Summary 

We have presented a unified theoretical description of the many-body states 
of a few electrons (holes) or few excitons confined in coupled quantum dots. 
In these systems, interdot coupling controls the competition between kinetic 
energy and Coulomb interactions, and can reach regimes far out of those ac- 
cessible in natural molecules. The resulting ground state is therefore very 
different for different values of the coupling: We have shown that a system 
of few electrons is characterized by different spin configurations depending on 
the inter-dot coupling, and we have extensively discussed the marked varia- 
tions arising in the electron-electron correlation functions. In the case of two 
electrons and two holes we have identified the ground state corresponding to 
both pairs localized in one of the dots (weak coupling) or distributed in both 
dots (strong coupling). Tuning such phases by external fields is possible, and is 
found to induce novel quantum effects that will be described elesewhere [72]. 

Manifestations of transitions between such phases in addition or optical spec- 
tra are expected to lead to a direct experimental verification of many-body- 
theory predictions, and to the experimental control of the many-body states 
in nanoscale devices. 
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M r eh r el)4>h(rhh r H)- 
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Fig. 1. Dependence on the inter-dot distance d of: (a) the kinetic (Ef.) 
and (b) Coulomb {V ee ) energies of three prototypical two-electron states 
|1) = \a g 1,a g I), |2) = (\a g ],a u |) - \a g [ , a u ]))/V2 (singlet states) and 
1 3) = {\a g j, o u I) + \a g I, a u |))/y / 2 (triplet state); the diamond symbol refers to the 
real ground state \tp). In panel (c) we plot the spatially averaged pair-correlation 
function p{z) = JfP (x, y, z; xq, yo, zq) dxdy corresponding to these states. The 
coordinates of the fixed particle (represented by the black circle) are xq = , yo = 
and zq = —7.5 nm ( z = is in the middle of the inter-dot barrier). The contin- 
uous grey line gives the profile of the z component of the confinement potential 
(barrier height Vq = 400 meV, in-plane confinement energy hujQ = 20 meV). Other 
parameters adopted are those typical of GaAs: m* = 0.067 m e , K r = 12.9. 
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Fig. 2. Angular pair correlation function P(g, <p, z; go,(fo, Zo) vs the azimuthal angle 
ip for N = 2: all other parameters g, go, z, zq are fixed, with g = go. Here go, zq, 
z correspond to the average value of the in-plane radius and the maxima along z 
of the single-particle density, respectively. We use m* = 0.067m e , n r =12.9, L = 10 
nm, Vq = 400 meV, d = 1 nm. The inset is a contour plot of P(g, z; go, zq) in the 
xy plane (in units of to), for the triplet excited state with Hujq = 20 meV. 
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Fig. 3. Ground state energy vs d for different number of electrons. Some excited 
states are also depicted, toghether with their term in the Molecular Spectroscopy 
notation. The bottom panel pictorially shows the single-particle configurations that 
have the largest weight in the three different many-body ground states for N = 4. 
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Fig. 4. Spin-resolved pair correlation function P^^(g, z; qq, zq) for N = 4 for three 
different values of d (columns). The position of a spin-up electron (qo,zq) is fixed 
in one dot as in Fig. 2 (black bullet), and the contour plots of the top (bottom) 
row, with z = zq {z ^ zq) fixed, correspond to the probability of measuring another 
spin- up electron on the same (other) dot in the xy plane (£q units). 
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Fig. 5. Plot of the functions p(z e ) (continuous line) and p(zh) (dashed line), com- 
puted for the excitonic ground states at d = 1 nm (a) and d = 3 nm (b). The defi- 
nition of p(z) and the positions of the fixed particles (either electron or hole) are as 
in Fig. 1. The columns in the insets represent the square moduli of the coefficients 
associated with the states |1) = \a g |,cr 9 j) (grey column) and |2) = \a u 1,cr u I) 
(white column): the excitonic ground states are given, to a good degree of approx- 
imation, by the superposition of these two states. At d = 1 nm |ci| 2 = 0.871 and 
| C2 1 2 = 0.066, while at d = 3 nm |ci| 2 = 0.577 and \c2\ 2 = 0.344; minor contributions 
arise from occupations of higher-energetic single-particle states. 
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Fig. 6. The two spin configurations (S e = 0,Sh = 0) and (S e = l,Sh = 0) are 
considered at d = 1 nm. The plots represent the spatial distribution of a carrier of 
spin orientation as specified by the subscript of the z-coordinate, given the position 
of another carrier is fixed (whose type and spin orientation are drawn close to the 
corresponding black circle). The position of the fixed particles is the one adopted 
in Fig. 1. 
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Fig. 8. Plots (a), (b) and (c) represent the d-dependence of the average values of 
the different contributions to the Coulomb energy for the eigenstates of lowest total 
energy associated to each of the four spin configurations (S e ,Sh)- The curves asso- 
ciated to the (1,0) and (1, 1) spin arrangements are identical in (a), and, together 
with (0, 1), in (c). In (d) we plot the total energies. 
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